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Abstract. The physics of non-zero temperature dynamics and transport 
near quantum-critical points is discussed by a detailed study of the O(N)- 
symmetric, relativistic, quantum field theory of a iV-component scalar field 
in d spatial dimensions. A great deal of insight is gained from a simple, exact 
solution of the long-time dynamics for the N = 1 d = 1 case: this model de- 
scribes the critical point of the Ising chain in a transverse field, and the dy- 
namics in all the distinct, limiting, physical regions of its finite temperature 
phase diagram is obtained. The N = 3, d = 1 model describes insulating, 
gapped, spin chain compounds: the exact, low temperature value of the spin 
diffusivity is computed, and compared with NMR experiments. The N = 3, 
d = 2, 3 models describe Heisenberg antiferromagnets with collinear Neel 
correlations, and experimental realizations of quantum-critical behavior in 
these systems are discussed. Finally, the N = 2, d = 2 model describes the 
superfluid-insulator transition in lattice boson systems: the frequency and 
temperature dependence of the the conductivity at the quantum-critical 
coupling is described and implications for experiments in two-dimensional 
thin films and inversion layers are noted. 
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1. Introduction 

Consider a quantum system on an infinite lattice described by the Hamilto- 
nian Tt(g), with g a dimensionless coupling constant. For any reasonable g, 
all observable properties of the ground state of TL will vary smoothly as g is 
varied. However, there may be special points, like g = g c , where there is a 
non-analyticity in some property of the ground state: we identify g c as the 
position of a quantum phase transition. In finite lattices, non-analyticities 
can only occur at level crossings; the possibilities in infinite systems are 
richer as avoided level crossings can become sharp in the thermodynamic 
limit. In this paper, I will restrict my discussion to second order quantum 
transitions, or transitions in which the length and time scales over which the 
degrees of freedom are correlated diverge as g approaches g c . As I will dis- 
cuss below, any such quantum transition can be used to define a continuum 
quantum field theory (CQFT): the CQFT has no intrinsic short-distance 
(or ultraviolet) cutoff. The main purpose of this paper is to discuss some 
properties of 7~C(g) at finite temperatures (T) in the vicinity of g = g c , in 
the context of some simple, but experimentally important, models. These 
studies are equivalent to a determination of the finite T crossovers of the 
associated CQFTs. 

We begin by stating some basic concepts on the relationship between 
quantum critical points and CQFT's [1, 2, 3]. As correlations become long 
range in time in the vicinity of the critical point, every system must be 
characterized by an experimentally measurable energy scale, A which van- 
ishes at g = g c . Convenient choices are an energy gap, if one exists, or a 
stiffness of an ordered phase to changes in the orientation of an order pa- 
rameter. In the models we shall consider here, A vanishes as a power-law 
as g approaches g c : 

A~A.\g-g c r (1) 

where A is an ultraviolet cutoff, measured in units of energy, z is the dy- 
namic exponent, and v is the correlation length exponent [4, 5]. From the 
perspective of a field theorist, the CQFT associated with the quantum crit- 
ical point is now defined by taking the limit 

bark — > oo at fixed A; from (1) we see that, because zv > 0, it is possible 
to take this limit by tuning the bare coupling g closer and closer to the 
critical point as A increases. (A condensed matter physicist would take the 
complementary, but equivalent, perspective of keeping A fixed but moving 
closer to criticality by lowering his probe frequency u ~ A). Assuming 
the A —> oo limit exits, the resulting CQFT then contains only the energy 
scale A. At non-zero temperatures, there is a second energy scale ksT; 
its thermodynamic properties will then be a universal function of the only 
dimensionless ratio available — A/ksT. 
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Figure 1. Schematic phase diagram as a function of the coupling constant of the quantum 
Hamiltonian g, and the temperature T. The quantum critical point is at T — 0, g = g c . 
The dashed lines indicate crossovers. There may finite temperature phase transitions in 
either of the two low temperature regimes. The symbol A represents a characteristic 
energy scale of the T = theory; it vanishes at g approaches g c according to (1). 

Our study will find two regimes with very different physical properties, 
as sketched in Fig 1: 

(i) The low temperature region fcgT< A 

There are actually two regimes of this type, one on either side of g c . Cor- 
relations in this region are similar to those of the T = ground state. The 
low temperature creates a small density of excitations which can sometimes 
have significant effects at very long scales. 

(ii) The high temperature region A <C fcgT 

As we are discussing universal properties of the CQFT, it is implicitly as- 
sumed that all energy scales, including ksT, are smaller than the upper 
cutoff A which has been sent to infinity; so we also require that fcsT< A. 
The thermal energy, fegT, sets the scale for all physical phenomena in this 
region, and the system behaves as if it's couplings are at the g = g c critical 
point. We shall devote much attention to the unfamiliar and unusual prop- 
erties of this region. It is perhaps worth noting explicitly why the high T 
limit of the CQFT can be non-trivial. A conventional high T expansions of 
the lattice model H proceeds with the series 

TVe -«/^ = m _^ w+ _i_ ltw2 + ... (2) 

The successive terms in this series are well-defined and finite because of 
the ultraviolet cutoffs provided by the lattice. Further, the series is well- 
behaved provided T is larger than all other energy scales; in particular we 
need k^T 3> A. In contrast, the CQFT was defined by the limit A — > oo 
at fixed ksT, A, and, as already stated, the high T limit of the CQFT 
corresponds to the intermediate temperature range A <C ksT <C A of 
the lattice model. It is not possible to access this temperature range by 
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an expansion as simple as (2), and more sophisticated techniques, to be 
discussed here, are necessary. 

All of our explicit computations will be with models described by the 
same CQFT. This is the relativistic CQFT in d spatial dimensions with 
imaginary time (r) action 



s = J dTjd d xt [ - [(d T <p a ) 2 + c 2 (vm 2 

+K c + ( 5 -5cM] + ^(^) 2 }.(3) 

Here (j) a is a real scalar field, the index a = 1 ... N is implicitly summed 
over, and the action has O(N) symmetry. This CQFT has a "Lorentz" in- 
variance with c the velocity of "light" , and as a result the dynamic critical 
exponent z = 1. The bare "mass" term has been written as + (g — g c ) 
so that the T = quantum critical point is at g = g c , and no measures the 
strength of the quartic non-linearity. The action S can also be interpreted 
as the Gibbs weight of a classical statistical mechanics problem in d + 1 di- 
mensions (with the r dimension of finite extent h/ksT); indeed it is nothing 
but the standard, thoroughly-studied 4> 4 theory which is the corner-stone 
of the well established theory of classical critical phenomena [1]. It might 
then appear that we can simply carry over these results to the case of the 
quantum-critical point, and our job is relatively straightforward. This is far 
from being the case. The fundamental reason is that all dynamic experimen- 
tal measurements are in real time (t), and we are particularly interested in 
the long time limit t S> U/ksT. While, in principle, information on these 
long-time correlations is related by analytic continuation to imaginary time 
correlations in the domain < r < H/ksT, in practice the continuation is 
an ill-posed problem, and essentially impossible to carry out. In particular, 
it has been shown [6, 7, 8] that the operations of expansion in e = 3 — d or 
1/N (which are the only non-numeric tools available for analyzing S), and 
of analytic continuation do not commute. It is essential that the theory of 
the dynamic and transport properties of S be formulated directly in real 
time, and here we shall review recent progress in this direction. It is sober- 
ing to note that there remain open questions on experimentally important 
observables even for the simple model S, and we shall also note them here. 

A central concept in our description of the dynamic properties of S is 
that of the phase relaxation time, r^. This is defined as the time over which 
the wave- functional of the CQFT retains phase memory. We shall find that 
in the regions of Fig 1 

t v ~ h/k B T in the "High T" region 

r v > h/k B T in the "Low T" regions (4) 
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These relations will be shown to apply to S, but are expected to be far more 
general. The missing constant in the first of these relations is a universal 
number which depends upon the precise definition of r^. All of the physical 
properties of S will show significant crossovers at times of order t v which 
we shall describe. 

2. The Ising chain in a transverse field 

We will begin by obtaining exact results for the dynamic scaling functions 
of S for the case d = 1, N = 1. The quantum physics is more transparent in 
a lattice Hamiltonian formulation, where the properties of S are expected 
to be equivalent to those of the Ising chain in a transverse field (this may 
be shown in a manner similar to Section 3.2 which considers N > 1). The 
degrees of freedom of the Ising model are spins af on the sites, i, of a chain. 
In addition to their usual exchange interaction J, there is a transverse field 
of strength gj which is responsible for the quantum dynamics; thus the 
Hamiltonian is 

F; = -J^(K + «i) (5) 

i 

Like S for N = 1, Hj possesses a global Zi symmetry: it is invariant under 
a unitary transformation performed by the operator \\ i af under which 

(6) 

Notice that this global Ising Z<i symmetry is present in the presence of 
the transverse field. A longitudinal field, coupling to a z would break the 
Zi symmetry. The action of this symmetry correctly indicates that the 
operator correspondence a z ~ 4> maps the long distance correlators of S 
and Hi. 

2.1. LIMITING CASES 

We begin by examining the spectrum of Hi under strong (g S> 1) and 
weak (<7 <C 1) coupling [9]. The analysis is relatively straightforward in 
these limits, and two very different physical pictures emerge. The exact 
solution, to be discussed later, shows that there is a critical point exactly 
at g = 1, but that the qualitative properties of the ground states for g > 1 
(g < 1) are very similar to those for g » 1 (g <C 1). One of the two 
limiting descriptions is therefore always appropriate, and only the critical 
point g = 1 has genuinely different properties. 

2.1.1. Strong coupling g 3> 1 

Let and | |)j denote the eigenstates of af. Then |±)j = (|T)j± I l)i)/V% 
are the eigenstates of af. Then at g = oo the ground state of Hi is clearly 



a- 



-a, 
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determined by the transverse field term to be 

io>=ni+>< w 

i 

The values of of on different sites are totally uncorrected in this state, and 
so (0|of <rj|0) = Sij. Perturbative corrections in l/g will build in correlations 
in a z which increase in range at each order in l/g; for g large enough these 
correlations are expected to remain short-ranged, and the a z correlator to 
decay exponentially with separation. There is thus no magnetic long-range 
order and this state is a "quantum paramagnet". Notice that this state is 
invariant under the Z2 symmetry described above. 

What about the excited states ? For g = 00 these can also be listed 
exactly. The lowest excited states are 

i<>= i->i ni+>i. ( g ) 

obtained by flipping the state on site i to the other eigenstate of a x . All 
such states are degenerate, and we will refer to them as the "single-particle" 
states. Similarly, the next degenerate manifold of states are the two-particle 
states \i,j), obtained by flipping the states at sites i and j, and so on to the 
general n-particle states. To leading order in l/g, we can neglect the mixing 
between states between different particle number, and just study how the 
degeneracy within each manifold is lifted. For the one-particle states, the 
exchange term in Hj leads only to the off-diagonal matrix element 

(i\Hj\i + l) = -J (9) 

which hops the 'particle' between nearest neighbor sites. As in the tight- 
binding models of solid state physics, the Hamiltonian is therefore diago- 
nalized by going to the momentum space basis 

l fc > = 7=E efcl l*> (io) 

where iV is the number of sites. This eigenstate has energy (we have choosen 
an overall constant in Hj to make the energy of the ground state zero) 

e k = Jg (2- (2/g) cos ka + 0(l/g 2 )) (11) 

where a is the lattice spacing. The lowest energy one-particle state is there- 
fore at Eq = 2g — 2 J 

Now consider the two-particle states. As long as the two particles are 
well separated from each other, the eigenstate is formed simply by taking 



7 



the tensor product of two single particle eigenstates. However these parti- 
cles will collide, which will be described by a 5 matrix. If we exclude the 
possibility of two-particle bound states (which do not occur here), the total 
energy of the state is determined by the configuration where the particles 
are well separated, and is simply the sum of the single particle energies. 
Thus the energy of a two-particle state with total momentum k is given by 
Ek = £ki + £ fc 2 w here k = k\ + ki- Notice that for a fixed k, there is still an 
arbitrariness in the single particle momenta k±^ and so the total energy 
can take a range of values. There is thus no definite energy momentum re- 
lation, and we have instead a 'two-particle continuum'. It should be clear, 
however, that the lowest energy two-particle state in the infinite system 
(its "threshold") is at 2e$. Similar considerations apply to the n-particle 
continua, which have thresholds at ne$. 

At next order in 1/g we have to account for the mixing between states 
with differing numbers of particles. Non-zero matrix elements like 



lead to a coupling between n and n + 2 particle states. It is clear that these 
will renormalize the one-particle energies St- However qualitative features 
of the spectrum will not change, and we will still have renormalized one- 
particle states with a definite energy-momentum relationship, and renor- 
malized n > 2 particle continua with thresholds at heq. 

The same expansion in powers of 1/g can also be used to compute the 
two-particle S matrix. If we consider the collision of two-particles with 
small momenta p and p' , then by conservation of energy and momentum 
there can only be two particles in the final state, and the momenta of these 
particles remains p and p'. An elementary calculation to order 1/g shows 
then that 



It can be shown that this result holds in limit of small p, p' to all orders in 
1/g for a large class of models like Hj with further neighbor exchange. The 
particles under consideration, being simple spin flips, are evidently bosons 
which experience a short-range repulsive potential. In d = 1 any weak 
repulsive potential appears arbitrarily strong in the limit of low velocities, 
and leads to the 'unitarity limit' phase shift of n, which is responsible for 
(13). For the particular nearest neighbor model Hj, it can be shown that 
there is no particle production in collisions at any incoming momenta p, p' , 
and that (13) holds to all orders in 1/g for all p, p'. This will be shown and 
exploited in Section 2.2. 

The spectrum described above has simple, but important, consequences 
for the dynamic spin susceptibility x(A;,iu; ra ). This is defined in imaginary 



(0|#/|m + 1) = - J 



(12) 




= -1. 



(13) 
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Figure 2. Schematic of the spectral density ~x'(k, w) of Hi as a function of ui at T = 
and a small k. There is a quasiparticle delta function at uo — et, and a three-particle 
continuum at higher frequencies. 

time as the Fourier transform into momentum and frequency of the a z 
correlator: 

1 /T 

x(fc,iw„)=jf dr J dx(a z (x, ir)a z (0,0)) e ^ kx ~^ (14) 

Its spectral density x"{Pi u ) is the imaginary part of the real frequency 
x(k,u), and it given by 

x"(k,co) = vr£ | <0|<7*(fc)|a> \ 2 8(lo - E a ) (15) 

a 

where the sum over a extends over all the eigenstates of Hj with energy E a . 
The eigenstates and energies described above allow us to simply deduce the 
qualitative form of x"iPi w ) which is sketched in Fig 2. The operator a z flips 
the state at a single site, and so the matrix element in (15) is non-zero for 
the single particle states: only the state with momentum p will contribute, 
and so there is an infinitely sharp delta function contribution to x"{k,uj) ~ 
5(u — £;-)■ This delta function is the "quasiparticle peak" and its co-efficient 
is the quasiparticle amplitude. At g = oo this quasiparticle peak is the entire 
spectral density, but for smaller g the quasiparticle amplitude decreases and 
the multiparticle states also contribute to the spectral density. The mixing 
between the one and three particle states discussed above, means that the 
next contribution to x"(p> a; ) occurs above the 3 particle threshold uj > 3eq; 
because there are a continuum of such states, their contribution is no longer 
a delta function, but a smooth function of omega (apart from a threshold 
singularity), as shown in Fig 2. Similarly there are continua above higher 
odd number particle thresholds; only states with odd numbers of particles 
contribute because the matrix element in (15) vanishes for even numbers 
of particles. 
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2.1.2. Weak coupling g <C 1 

Now the energy is dominated by the exchange term. There are two de- 
generate ground states at g = with the spins either all up or down (in 
eigenstates of a z ): 

it> = n i ?>* u> = n u>i ( ie ) 

i i 

Turning on a small g will mix in a small fraction of spins of the opposite 
orientation, but the degeneracy will survive as the two states are related to 
each other by the global Z% symmetry noted above (6). A thermodynamic 
system will always choose one or the other of the states as its ground states 
(which may be preferred by some infinitesimal external perturbation), and 
hence the Z2 symmetry will be spontaneously broken. The correlations of 
the magnetization a z have an infinite range in either state as 

lim K(x,0K(0,0)) = iV V0 (17) 

|x|— >oo 

The quantity Nq is the spontaneous magnetization, and equals (a z ) in ei- 
ther of the two ground states. All of the statements made in this paragraph 
clearly hold for g = 0, and will hold for some g > provided the perturba- 
tion theory in g has a non-zero radius of convergence. The exact solution 
of the model to be discussed later will verify that this is indeed the case. 

The excited states can be described in terms of an elementary domain 
wall (or kink) excitation. For instance the state 

" " " lt)i lt)i+l ll)i+2 ll)i+3 U)i+4 IT)i+5 IT)i+6 ' ' ' 

has domain walls, or nearest neighbor pairs of antiparallel spins, between 
sites i + i + 2 and sites i + 4, i + 5. At g = the energy of such a state is 
clearly 2Jxnumber of domain walls. The consequences of a small non-zero 
g are now very similar to those due to \ j g corrections in the complementary 
large g limit: the domain walls become "particles" which can hop and form 
momentum eigenstates with excitation energy 

e k = j(2-2gco S (ka) + 0(g 2 )) . (18) 

The spectrum can be interpreted in terms n-particle scattering states, al- 
though it must be emphasized that the interpretation of the particle is now 
very different from that in the large g limit. Again, the perturbation theory 
in g only mixes states which differ by even numbers of particles, although 
now the matrix element in (15) is non-zero only for states a with an even 
number of particles; these assertions can easily be checked to hold in a per- 
turbation theory in g. So x"(Pj u ) wm now have a pole at p = 0, uj = + , 
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from the term in (15) where a = one of the ground states, indicating the 
presence of long-range order. Further, there is now no single particle con- 
tribution, and the first finite u spectral density is the continuum above the 
two particle threshold. The absence of a single particle delta function in 
this case is a very special feature of the d = 1, N = 1 model, and is not 
expected to hold in higher d. 

The S matrix for the collision of two domain walls can now be computed 
in a perturbation theory in g, and as in the strong-coupling 1/g expansion, 
we find that there is no particle production, and S pp i = — 1 to all orders in 
9- 

2.2. EXACT SPECTRUM AND CONTINUUM THEORY 

The qualitative considerations of the previous section are quite useful in 
developing an intuitive physical picture. We will now take a different route, 
and set up a formalism that will eventually lead to an exact determination 
of many physical correlators; these results will vindicate the approximate 
methods for g > 1, g < 1, and also provide an understanding of the novel 
physics at g = 1. 

The central idea is the application of the Jordan- Wigner transforma- 
tion [10, 11]. We map the |+), |— ) states on each site to the Fock space of 
spinless fermions which can have occupation numbers 0, 1 on each site. The 
operator representation of the mapping is 

of = l-2e\* 

a* = (-l/2)n(l-2ctc,-)(ci + 4) (19) 

j<i 

where q is the annihilation operator for the fermion at site i. Inserting (19) 
into Hj we find that the resulting Hamiltonian is quadratic in the fermions 

Hi = - j J2 (4 c m + 4+i c i + 44+i + c *+iCj - 2gc\ci - gj . (20) 

i 

This Hamiltonian can be diagonalized in momentum space by a Bogoliubov 
transformation. The fermionic quasiparticle operators are 7& (—it /a < k < 
ir /a) and the diagonal form of Hj is 

Hi = Y, £ k{ibk-im (21) 

where 

e k = 2 J (l + g 2 - 2g cos(fca)) 1/2 (22) 
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is the single particle energy. As > 0, the ground state, |0), of Hi has 
has no 7 fermions and therefore satisfies 7fc|0) = for all k. The excited 
states are created by occupying the single particle states; they can clearly 
be classified by the total number of occupied states, and a n-particle state 
has the from 7^7^ • • "7^ |0), with all the ki distinct. 

The above structure of the spectrum confirms the approximate consid- 
erations of Section 2.1. We have now found that the particles are in fact free 
fermions, and two fermions will not scatter even when they are close to each 
other; alternatively they can be considered as hard-core bosons which have 
an S matrix which does not allow particle production, and which equals 
— 1 at all momenta. We shall find it much more useful to take the latter 
point of view, as the bosonic particles have a simple, local, interpretation 
in terms of the underlying spin excitations: for g S> 1 the bosons are simply 
spins oriented in the |— ) direction, while for g <C 1 they are domain walls 
between the two ground states. The fermionic representation is useful for 
certain technical manipulation, but the bosonic point of view is much more 
useful for making physical arguments, as we shall see below. 

The excitation energy in (22) is non-zero and positive for all k pro- 
vided g 7^ 1. The energy gap, or the minimum excitation energy is always 
at k = 0, and equals 2J|1 — g\. This gap vanishes at g = 1, and it is nat- 
ural to expect that g = g c = 1 is the phase boundary between the two 
qualitatively different phases discussed in Section 2.1. Precisely at g = 1, 
fermions with low momenta can carry arbitrarily low energy, and therefore 
must dominate the low temperature properties. These properties suggest 
that the state at g = 1 is critical, and there is a CQFT which describes 
the critical properties in its vicinity. As the important excitations are near 
k = 0, we define the continuum Fermi field 

*(xi) = -±=* (23) 
V a 

We express i7j in terms of and expand in spatial gradients. This gives 
the fermionic Lagrangean of the CQFT believed to be equivalent to S for 
d= 1, N= 1: 

The field * is now implicitly assumed to be a function of space and imagi- 
nary time (r). The coupling constants in (24) are given for Hi by 

A = 2.7(1 -0) c = 2Ja. (25) 

However, the relations (25) are very specific to the solvable model 77/. For 
more complicated, non-solvable, Ising models which have a similar naive 
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continuum limit (e.g. models with second neighbor exchange), the values of 
A and c appearing in the continuum quantum field theory cannot be deter- 
mined exactly. Indeed A and c are parameters which depend upon details 
of the microscopic Hamiltonian, i.e. they have a non-universal dependence 
upon a coupling constant like g, as can be checked by a simple estimate of 
perturbative fluctuation corrections due to irrelevant operators. 

Comparing the dependence of A on g in (25) with (1) we conclude that 
the exponent zv = 1. Further at the critical point A = 0, the excitation 
energy ~ k, which fixes the dynamic critical exponent z = 1. 

The continuum theory Cj can be diagonalized much like the lattice 
model H[, and the excitation energy now takes a "relativistic" form 

e k = (A 2 + c 2 k 2 ) 1/2 (26) 

which shows that |A| is the T = energy gap (we will choose the sign 
of A to be different on the two sides of the critical value of g), and c is 
the velocity of the excitations, both measurable quantities. The form of 
correctly suggests that Lj is invariant under Lorentz transformations. This 
can be made explicit by writing the complex Grassman field ^ in terms of 
two real Grassman fields, when the action becomes what is known as the 
field theory of Major ana fermions of mass A/c 2 [12]; we will not explicitly 
display this here. 

While the continuum field ^ has no direct physical interpretation, its 
correlators are easily computed at all g, including the critical point g = 1. 
To gain some insight about this point, we explicitly display a certain A = 
two-point correlator for T > 0. We find for imaginary time r > 

»(x,*).t(0,0» = \ Jl (e*<>/-> 

+ ',,_,A (27. 



4c/ Vsin(7rT(r — ix/c)) sin(7rT(r + ix/c)) 

We are now using units in which h = ks = 1, and will continue to do so in 
the remainder of the paper. At T = 0, (27) simplifies to 

*(x,r)*t( ,0)) = -!- (—^— + —^—) , (28) 

' 47T V CT — IX CT + IX J 

when we notice that the transformation 



CT 



c firT \ 

± ix — > — — sin (cr ± ix) (29) 

ttT \ c ) 



connects the T = and T > results. This transformation is actually a 
very general property of the critical point, and connects all T = and 
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Figure 3. Finite T phase diagram of the d = 1 quantum Ising model, Hj, as a function of 
the coupling g and temperature T. There is a quantum phase transition at T — g = g c 
with exponents z = 1, v = 1. Long-range order (No = (cr z ) 7^ 0) is present only for T = 
and g < g c - There is an energy gap above the ground state for all g ^ g c . We use an 
energy scale A ~ g c — g such that the energy gap is |A|. The dashed lines are crossovers 
at |A| ~ T. A very similar phase diagram also applies to the d = 2 0(3) quantum rotor 
model discussed in Section 3.4. In this case there is no energy gap for g < g c , and the 
spin-stiffness p 3 is used as the energy scale to characterize the ground state. There is an 
energy gap A for g > g c . These energies vanish as p s ~ (g c — g) zv , A ~ (g ~ g € ) zv with 
z = 1, v » 0.7, and so the crossover phase boundaries at T ~ |ff"5ci z " will not be linear. 



T > correlators; it is a consequence of the conformal invariance [13] of Cj 
at A = 0. We will use (29) in an important way later. 

Correlators of a x can be constructed out of those of simple bilinears 
of the fermion operators, and we will not display them explicitly. More 
interesting, however, are the correlators of the order parameter a z ~ (p. 
Computing just the equal time two-point correlator, or even simply the 
value of its scaling dimension, dim[cr 2 ], involves a rather lengthy and in- 
volved computation, which will not be discussed here. Rather, in the next 
section, we will quote a recently obtained technical result, and then proceed 
to obtain dynamic correlators of a z by simple physical arguments. 

2.3. FINITE TEMPERATURE CROSSOVERS 

The purpose of this section is to describe the dynamics of the order param- 
eter in the different regions of the T > phase diagram sketched in Fig 3; 
this diagram follows from the considerations in Section 1, the expression 
for A in (25), and computations to be described below. As noted earlier, 
here we will quote just one technical result on the equal-time two-point 
correlator of a z for T > which was obtained [14, 15, 16] using the fermion 
mapping. We will then show, following the recent work of Ref. [17], that 
unequal-time, T > 0, correlations can be obtained by simple physical ar- 
guments that rely on the bosonic picture of the excitations developed in 
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Section 2.1 using the large and small g expansions. It is worth noting that 
very sophisticated methods [18, 19] relying on the fermion mapping have 
not so far succeeded in obtaining any explicit results for unequal time corre- 
lations for T > (there are some results [20] on time-dependent correlators 
at T = oo which are non-universal and unrelated to the CQFT of interest 
here) , 

The equal-time, long-distance result we need is [16] 

(a z (\x\ -► oo,t = 0)<7*(0,0)) = ZT^GKA/Tjexp (-^-F^A/T) 

(30) 

where t is real time, Z is non-universal constant, and Fj(s) and Gi(s) are 
universal scaling functions. The crucial property of (30) is the prefactor 
of T 1 / 4 . As T is an energy which scales as inverse time, and the dynamic 
exponent z = 1, this allows us to conclude that the scaling dimension of a z 
(0) is 

dim[a z ] = 1/8. (31) 
From (30), we can also define the correlation length £ which obeys 

c'^ Fl (*) (32) 



c \T, 

The exact, self-contained expression for the universal function Fj is [16] 

F I (s) = - cfo/lncoth ^ ^ ' +\s\0(-s). (33) 
vr Jo 2 

The s > (s < 0) portion of F/ describes the ordered (disordered) side. 
Despite appearances, the function Fj(s) is smooth as a function of s for all 
real s, and is analytic at s = 0. The analyticity at s = is required by the 
absence of any thermodynamic singularity at finite T for A = 0. This is 
a key property, which was in fact used to obtain the answer in (33). The 
exact expression for the function Gi(s) is also known 

ta0 , w = f d JL \( d JM) 2 _ i] + r * f^M) 2 , (34) 

J s y [\ dy J 4J 7i y \ dy J 

and its analyticity at s = follows from that of Fj. For the solvable model 
Hi, we chose the overall normalization of Gj such that Z = J^ 1 / 4 . In 
general, the value of Z is set by relating it to an observable, as we will 
show below. Also note that Z has no dependence on A, and is therefore 
non-singular at the quantum critical point. 
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Armed with the above knowledge, we can write down the full scaling 
form for the time-dependent o z correlator, which applies to the lattice 
model in the limits A ~ J ^ oo, a^Oat fixed A, c and T 

<<7*(MK(0,0)> = ZT 1 ^ (^,7V,^) (35) 

where $j is a universal function which is analytic as a function of its third 
argument s = A/T on the real s axis. The result (30) obviously specifies 
<&/ for large Tx/c and t = 0. 

The following subsections will describe the unequal-time form of in 
the limiting regions of Fig 3: they are associated with the limits s — ► oo 
(renormalized classical), s — > — oo (quantum disordered), and s = (high 
T) of Fj, Gj which will also be noted below. 

2.3.1. Low T on the ordered side, A > 0, T <C A 

This is the "renormalized classical" [21] region of Fig 3, and the reasons for 
this name will become clear below. 

Assuming that it is valid to interchange the limits T — > and x — > oo 
in (30), we can use the limiting values Fj (oo) = 0, Gj(s — > oo) = s 1 / 4 to 
deduce that (recall (17)): 

A 2 = lim (a z (x,0)a z (0,0)} = ZA 1/4 at T = 0. (36) 

\x\— »oo 

Thus, as claimed earlier, there is long-range order in the g < 1 ground state 
of Hi, with the order parameter Nq = (a z ) = Z 1 / 2 A 1 / 8 (this relates the 
value of Z to a physical observable). For small T < A, we obtain from the 
large s behavior of Fi(s) (see (33)) that 

(<T z (x,0)<j*(p,0))=N$e-Wt° large \x\, (37) 
where the correlation length 

{ -(^)"V- 

is finite at any non-zero T, showing that long-range order is present only 
precisely at T = 0. We have put a subscript c on the correlation length 
to emphasize that the system is expected to behave classically in this low 
temperature region [16]. The excitations above the ground states consists of 
particles (the kinks and anti-kinks of Section 2.1) whose mean separation (~ 
£ c ) is much larger than their de Broglie wavelengths (~ (c 2 / AT) 1 / 2 , as the 
mass of these particles = A/c 2 ), which is precisely the canonical condition 
for the applicability of classical physics. It is also reassuring to note that 
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(37) is precisely the form of equal-time correlations in the classical Ising 
model at low T. The prefactor Nq is the true ground state magnetization 
including the effects of quantum fluctuations, and this is the reason for the 
adjective "renormalized" in the name for this region. 

We now show, following Ref [17], that it is possible to give a simple 
physical interpretation for the precise value of £ c in (38). The energy of 
a domain wall with a small momentum k is A + c 2 k 2 /2A, and therefore 
classical Boltzmann statistics tells us that their density, p, is 

P = /f^ (a+AV2AVT =(H) 1/2 ^ i/T <») 

Comparing with (38), we see that £ c = 1/2/?. We now show that this result 
follows immediately if we assume that the domain walls are classical point 
particles, which are distributed independently with a density p. Consider 
a system of size L » \x\, and let it contain M thermally excited particles; 
then p = M/L. Let q be the probability that a given particle will be between 
and x. Clearly, q = \x\/L. The probability that a given set of j particles are 
the only ones between and x is then q^(l — q) M ~^: as each particle reverses 
the orientation of the ground state, in this case a z (x, 0)a z (0, 0) = Nq (— l)- 7 . 
Summing over all possibilities we have 

M M' 
(a z (x,0)a z (0,0)) = N 2 g(-l) V(l - ^ j[{M '_ j){ 

= N 2 (l-2q) N ^N 2 e~ 2qM = N 2 e~ 2 ^, (40) 

thus establishing the desired result. 

This classical picture can also be extended to compute unequal time cor- 
relations [17]. Let us explicitly consider the form of the correlation function 
needed 

(a z (x,t)a z (0,0)) = Tr (e~ HllT e iHlt a z e- iHlt <y^j /Z, (41) 

where Z = Tre~ H '/ T . We evaluate it in the double-time ('Keldysh') path 
integral formalism [22]. The path integral is over a set of trajectories moving 
forward in time, representing the operator e~ lHlt , and a second set moving 
backwards in time, corresponding to the action of e lHlt . In the semiclassical 
limit, stationary phase is achieved when the backward paths are simply the 
time reverse of the forward paths, and both sets are the classical trajecto- 
ries. An example of a set of paths is shown in Fig 4. Now observe that (i) 
the classical trajectories remain straight lines across collisions because the 
momenta before and after the collision are the same in one dimension; (ii) 
for each collision, the amplitude for the path acquires a phase S pp i along the 
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1 ► X 

Figure 4- A typical semiclassical contribution to the double time path integral for 
{cr z (x, t)a z (0, 0)). Full lines are thermally excited particles which propagate forward and 
backward in time. The ± signs are significant only for g < 1 and denote the orientation 
of the order parameter. For g > 1, the dashed line is a particle propagating only forward 
in time from (0,0) to (x,t). 



forward path and its complex conjugate along the backward path: the net 
factor for the collision is therefore (S'pp'l 2 = 1- These two facts imply that 
the trajectories are simply independently distributed straight lines, placed 
with a uniform density p along the x axis, with an inverse slope 

and with their momenta chosen with the Boltzmann probability density 

e -e k /T j p ( Fig 4 )_ 

Computing dynamic correlators is now an exercise in classical prob- 
abilities. The value of a z (0,0)a z (x,t) is the square of the magnetization 
renormalized by quantum fluctuations (Nq), times (—1) if the number of 
trajectories intersecting the dashed line in Fig. 4 is odd. This precisely the 
computation carried out above for the equal time correlations, but now it 
is easy to see from Fig 4 that for particles with velocity Vk 

(43) 



L 

Averaging over velocities, and then proceeding as in (40) we obtain one of 
our central results [17, 23] 

(a z (x,t)a z (0,0)) =iYo 2 exp(-| ^- e ~^ T \x - v k t\\ . (44) 

Notice that the equal-time correlator is Nq e~^^ c as before, and the equal- 
space correlator decays exponentially in time as N 2 e~\ t \/ Tc , but the general 
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behavior is more complicated. The classical relaxation time r c is a central 
quantity; remarkably, we find from (44) that the correlation time, r c , is 
independent of the functional form of and depends only on the gap: 

Wh e A ^ T , (45) 



2k B T 



where we have momentarily inserted the fundamental constants h, Jzb to 
emphasize the universality of the prefactor. As r c is also the time over which 
the phase coherence of the ground state is lost, we may identify the phase 
relaxation time of Section 1 as r v = r c . 

In the limit T<Awe are now able to completely specify the form of 
the scaling function $/ in (35). We find that <&/ collapses into a reduced 
scaling function <£/ c which is determined by classical thermal physics. In 
particular we have 

(a z (x, t)a z (0, 0)> = Nfa Ic (|, ^) , (46) 

with 

ln$ Jc (x,t) = -x erf (jj=) ~ te^l^. (47) 

Notice that the characteristic time r c and length £ c both diverge as ~ 
e A//T , and so we can define an effective classical dynamic exponent z c = 1 
(there is no fundamental reason why z c and z should have the same value) . 
Classical scaling forms like (46) have been discussed earlier [16], but it 
was incorrectly conjectured that the scaling functions would be those of 
the Glauber model [24] (which has z c = 2); Glauber dynamics does not 
conserve total energy and momentum, and these conservation laws have 
played a crucial role in the kinematic constraints on particle collisions. 

All of the results above have been compared with exact numerical com- 
putations and the agreement is essentially perfect [17], giving us confidence 
in the physical approach used to understand dynamic properties at T > 0. 

2.3.2. Low T on the disordered side, A < 0, T <C |A| 
This is the "quantum disordered" region of Fig 3. 

Now we need to take the s — > — oo limit of the functions Fj(s), Gi(s); 
from these limits we find 

ZT 

(cj 2 (x,0)cj 2 (0,0)) = 7^7§7I eHr|/f > M ->oo at fixed 0<T<|A|, (48) 
with the correlation length £ given by 



c \ 7TC 
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So correlations decay exponentially on a scale ~ c/|A|, and there is no 
long-range order. The equal time correlations at T = behave in a similar 
manner, although the limits T — > and \x\ — > oo do not commute for 
the prefactor of the exponential decay. The T = result is determined by 
the form-factor expansion technique [25], which we will not discuss here; it 
yields [26, 19] 

1 /2 

<^(x,0V 2 (0,0)> = ZIAI 1 / 4 ( ^p ) e- A N/ c , |x| - oo at T = 0. 

(50) 

The form factor expansion also yields the T = dynamic susceptibility 
(defined in (14)). The leading term is in fact precisely the quasiparticle 
pole at energy = (c 2 k 2 + A 2 ) 1 / 2 that was argued to exist in this phase 
in Section 2.1. We have 

2cZ\A\ 1 / 4: 

* k > u) = <W + A*-\u + iW + - ' T = ° (51) 

where S is a positive infinitesimal. The imaginary part of this gives the 
delta function sketched in Fig 2; the continuum of excitations above the 
three particle threshold come from higher order terms in the form factor 
expansion, and are represented by the ellipsis in (51). It can now be checked 
that the Fourier transform of (51) yields the leading term (50) in the equal 
time correlation function. 

The result (51) shows that the quasiparticle residue is 2ZA 1//4 (this 
is another relation between Z and a physical measurable, and along with 
(36) , it implies a relationship between the values of the residue and No as 
we approach the critical point from either side). The residue vanishes at 
the critical point A = 0, where the quasiparticle picture breaks down, and 
we will have a completely different structure of excitations. 

The above is an essentially complete description of the correlations and 
excitations of the quantum paramagnetic ground state. We now turn to the 
dynamic properties at T > 0. At nonzero T, there will be a small density of 
quasi-particle excitations which will behave classically for the same reasons 
as in Section 2.3.1: their mean spacing is much larger than their de Broglie 
wavelength. The collisions of these thermally excited quasi-particles will 
lead to a broadening of the delta function pole in (51): the form of this 
broadening can be computed exactly in the limit T <C |A| using a semi- 
classical approach similar to that employed for the ordered side [17]. The 
argument again employs a semiclassical path-integral approach to evaluat- 
ing the correlator in (41). The key observation now is that we may consider 
the operator a z to be given by 

a z (x, t) = 2Z\A\ 1 / 4 (4>(x, t) + ^(x, t)) + ... (52) 
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where tjr is the operator which creates a single particle excitation from the 
ground state, and the ellipsis represent multi-particle creation/annihilation 
terms which are sub dominant in the long time limit. This representation 
may also be understood from the g>l picture discussed earlier, in which 
the single-particle excitations where |— ) spins: the o z operator flips spins 
between the ±x directions, and therefore creates and annihilates quasipar- 
ticles. 

The computation of the nonzero T relaxation is best done in real space 
and time, so let us first write down the T = correlations in this represen- 
tation. We define K(x, t) to be T = correlator of the order parameter: 

K{x,t) = (a z (x,t)a z (0,0)} T =o 



I 



dk cZ\A\^ j kx _ Ekt 



2vr e fc 

K (\A\(x 2 -c 2 ty^ /c) m 



Z I A ! 1/4 v i \ 1^2 2.2x1/2 



7T 



where Kq is the modified Bessel function. This result is just the Fourier 
transform of (51. Note that for t > \x\/c, the Bessel function has imaginary 
argument and is therefore complex and oscillatory. 

Now we consider the evaluation of (41) in the semiclassical path- 

integral approach [17]. A typical set of paths contributing to the Keldysh 
path integral is still given by Fig 4, but its physical interpretation is now 
very different. The dashed line now represents the trajectory of a particle 
created at (0, 0) and annihilated at (x, t), and ± signs in the domains should 
be ignored. In the absence of any other particles this dashed line would 
contribute K(x,t) to (a z (x, t)a z (0, 0)). The scattering off the background 
particles (the full lines in Fig 4) introduces factors of the S'-matrix element 
Sppi; as the dashed line only propagates forward in time, the 5-matrix 
elements for collisions between the dashed and full lines (and only these) are 
not neutralized by a complex conjugate partner. Using the low momentum 
value Sppi = —1, we see that the contribution to (cr z (x, t)a z (0, 0)) equals 
(— l) ni K(x, t) where rig is the number of full lines intersecting the dashed 
line. The (— l) ne is precisely the term that appeared in the renormalized- 
classical, although for very different reasons. We can carry out the averaging 
over all trajectories as before, and obtain one of our main results [17] 

(a z {x, t)a z (0, 0)) = K(x, t) exp (- J ^- e - £ ^ T \x - v k t\j , (54) 

where K(x,t) is given by Eq. (53). 

The result (54) clearly displays the separation in scales at which quan- 
tum and thermal effects act. Quantum fluctuations determine the oscilla- 
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Figure 5. Theoretical and numerical results from Ref. [17] for the correlator 
{a z (x,t)a z (0,0)} of flj at x = 30 with J = 1, g = 1.1 (therefore A = -0.2), T = 0.1 
and so the system is in the quantum-disordered region of Fig 3. The numerical data was 
obtained for a lattice size L — 512 with free boundary conditions; it has a "ringing" at 
high frequency due to the lattice cutoff. The theoretical prediction is from the continuum 
theory prediction in Eq. (54) and is represented by the smoother curve. The envelope of 
the numerical curve fits the theoretical prediction well. 



tory, complex function K(x,t), which gives the T = value of the cor- 
relator. Exponential relaxation of spin correlations occurs at longer scales 
~ £ c , r c , and is controlled by the classical motion of particles and a purely 
real relaxation function. This relaxation leads to a broadening of the quasi- 
particle pole with widths of order £~ , r^ 1 in momentum and energy space. 
Note that the relaxation function is identical to that obtained in the renor- 
malized classical region, and so the expressions for £ c and r c are the same 
as before, upto the mapping A — * |A|. We can consider the presence of a 
quasi-particle pole in the response function as a consequence of quantum 
coherence in the ground state, and so the phase relaxation time, t^, over 
which this coherence is lost may be identified with r c ; we have therefore 



2k B T' 



JA|/fe s T 



(55) 



In Fig. 5 we compare the predictions of Eq. (54) with numerical results 
on a lattice of size L = 512. The theoretical curve was determined from the 
continuum expression for K(x,t), but the full lattice form for was used. 
The theory agrees well with the numerics; some differences are visible for 
small x, outside the light cone, but this is outside the domain of validity of 
(54). 
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2.3.3. High T, T > |A| 

Right at the critical point, A = 0, this regime extends all the way down 
to T = 0. We begin by writing the T = equal-time correlator of the 
continuum theory; from the scaling dimension of a z in (31), this must have 
the form 

(<r*(s, OK (0,0))~ 1 /4 atT = 0, A = 0, (56) 

We will now fix the prefactor in (56) using our earlier results. The key 
ingredient is our knowledge that the underlying continuum model Ci is 
is relativistically and conformally invariant at the T = critical point. 
Further we assume that the mapping (29) between correlators at T / 
from those at T = holds also for the two point correlator of a z Then from 
(56) we must have at T / 0, but A = 

K(x, irK(0, 0)) ~ T 1 / 4 m . (57) 

[sin(-7rT(r — ix/c)) sin(7rT(r + ix/c))] 1 

Notice that we are now working in imaginary time, as it is slightly more 
convenient for our purposes; the real time result can in this case be obtain by 
analytic continuation, as we have the exact functional form for < r < 1/T. 
Let us now use this result in the equal-time case in the regime xT/c S> 1. 
Precise results for this regime where quoted earlier in (30), where using the 
values Fj(0) = vr/4 (from evaluation of (33)) and Gj(0) = 0.858714569 . . . 
we have 

(a z (\x\ -► oo,r = OK (0,0)) = ZT 1/4 G/(0) exp (-^j^j at A = 

(58) 

Finally, comparing with (57) we obtain, for A = 

K(z, irK (0,0)> = ZT 1 / 4 — . (59) 

[sin(7rT(r — ix/c)) sin(7rT(r + ix/c))] ' 

As expected, this result is of the scaling form (35), and indeed completely 
determines the function for the case where its last argument is zero. 

Now let us turn to a physical interpretation of the main result (59). 
Consider first the case T = 0. By a Fourier transformation of the T = 
limit of (59) we obtain the dynamic susceptibility 

, r = o,A = o 

(60) 
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Figure 6. Spectral density, Imx(k,u))/u)Z, of Hi at its critical point g = 1 (A = 0) at 
T = 0, as a function of frequency uj, for a set of values of k. 



with 5 a positive infinitesimal. We plot Imx(k,uj)/uj in Fig 6. Notice 
that there are no delta functions in the spectral density like there were 
in the quantum disordered phase (Fig 2), indicating the absence of any 
well-defined quasiparticles. Instead, we have a critical continuum of excita- 
tions. However, the presence of sharp thresholds and singularities indicates 
that there is still perfect phase coherence, as there must be in the ground 
state. 

Now let us turn to non-zero T. We Fourier transform (59) to obtain 
x(k, iuj n ) at the Matsubara frequencies uj n and then analytically continue to 
real frequencies (there are some interesting subtleties in the Fourier trans- 
form to x(k,iu) n ) and its analytic structure in the complex co plane, which 
are discussed elsewhere[27]). This gives us the leading result for x(k,uj) in 
the high T region 



x(k,uj) = 



Zc G/(0)r(7/8) r 


f 1 ,uj + ck^ 
\16 +i 4vrT ) 


r 


f 1 . u — ck \ 
U6 +l 4vrT ) 


T 7 /4 4vr r(l/8) r 


(Vb .LO + ck\ 
V 16 + % 4vrT ) 


r 


(lh .uj — ck\ 
V 16 + 1 4vrT ) 



We show a plot of Imx/uj in Fig 7. This result is the finite T version of 
Fig 6. Notice that the sharp features of Fig 6 have been smoothed out on the 
scale T, and there is non-zero absorption at all frequencies. For uj,k T 
there is a well-defined peak in Imx/uj (Fig 7) rather like the T = critical 
behavior of Fig 6. However, for uj, k <C T we cross-over to the quantum 
relaxational regime [7] and the spectral density Imx/uj is similar to (but 
not identical) a Lorentzian around u = 0. This relaxational behavior can 
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Figure 7. The same observable as in Fig 6, T 11//4 ImY(fc, u>)/u>Z, but for This is 

the leading result for Imy for T 3> |A| i.e. in the high T region of Fig 3. All quantities 
are scaled appropriately with powers of T, and the absolute numerical values of both 
axes are meaningful. 



be characterized by a relaxation rate Tf> defined as 



r^ 1 = -i 



dlnx(0,u;) 



(62) 



u)=0 



(this is motivated by the phenomenological relaxational form x(0,u) = 
Xo/(l — zw/Tr + 0(u> 2 ))). The frequency scale on which the sharp features 
of Fig 6 have been smoothed out is also Tr, and so we may now identify 
the phase coherence time = 1/Tr. From (61) we determine: 

where we have returned to physical units. At the scale of the characteris- 
tic rate Tr, the dynamics of the system involves intrinsic quantum effects 
(responsible for the non-Lorentzian lineshape) which cannot be neglected; 
description by an effective classical model (as was appropriate in both the 
renormalized classical and quantum disordered regions of Fig 3) would re- 
quire that r^j <C hsT/fi, which is thus not satisfied in the high T region of 
Fig 3. Alternatively stated, the mean spacing between the thermally excited 
particles is now of order their de-Broglie wavelength, and so the classical 
thermal and quantum fluctuations must be treated at an equal footing. 

The ease with which our expressions for in (45,55,63) have been ob- 
tained bely their remarkable nature. Notice that we are working in a closed 
Hamiltonian system, evolving unitarily in time with the operator e _lff /*/^ ) 
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Figure 8. Crossovers as a function of frequency for the Ising model in the different 
regimes of Fig 3. The correlations in the two classical relaxational regimes are quite dif- 
ferent from each other. The "ordered" regime is in quotes, because there is no long-range 
order, and the system only appears ordered between spatial scales c/A and £ c . In the 
renormalized-classical and quantum-disordered regions \/t v ~ Te' A ' / ' T . 



from an initial density matrix given by the Gibbs ensemble at a temper- 
ature T. Yet, we have obtained relaxational behavior at low frequencies, 
and determined exact values for a dissipation constant. In contrast, in the 
theory of dynamics near classical critical points [28] , the relaxation dynam- 
ics is simply postulated in a rather ad hoc manner, and the relaxational 
constants are treated as phenomenological parameters to be determined by 
comparison with experiments. 

2.3.4. Summary 

The main features of the finite temperature physics of the quantum Ising 
model are summarized in Figs 3 and Fig 8. At short enough times or 
distances in all three regions of Fig 3, the systems displays critical fluc- 
tuations characterized by the dynamic susceptibility (60). The regions are 
distinguished by their behaviors at the low frequencies and momenta. In 
both the low T regimes of Fig 3 (renormalized classical and quantum disor- 
dered), the long time dynamics is relaxational and is described by effective 
classical models. The relaxation time, or equivalently, the phase coherence 
time, is of order (h/ kBT)e^ ener ^ g a P)/ fc s T ; anc [ j s therefore much longer 
than h/hsT. In contrast, the dynamics in the high T region is also relax- 
ational, but involves quantum effects in an essential way, as was described 
above. In this region the phase relaxation time is now of order h/ksT. 

Finally a few remarks about experiments. While there are no experimen- 
tal studies of the Ising chain in a transverse field, the quantum relaxational 
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dynamics of its critical point g = 1 is very closely related to the finite 
temperature dynamics of S = 1/2 Heisenberg spin chains. The spin corre- 
lators of the latter system are given by expressions very similar to those in 
Section 2.3.3, and the relaxational dynamics can has been experimentally 
probed in NMR experiments. There is good agreement between theory and 
experiments, and the reader is referred to some recent papers for further 
details [29, 30, 31]. 

3. 0(3) quantum rotors between one and three dimensions 

This section will examine S for N = 3 in d = 1, 2, 3: the N = 3 model 
describes the low energy physics of quantum antiferromagnets, and we will 
present a number of experimental applications. 

Many aspects of the physics will be similar to those described in Sec- 
tion 2 for the N = 1, d = 1 case. However there are some genuinely new 
features which shall be the focus of our attention. For N > 2, S has a contin- 
uous O(N) symmetry, and associated conserved charges, whose transport 
properties will be of interest to us. Further, for d > 2, there is a thermody- 
namic phase transition at a non-zero temperature: we shall briefly consider 
the subtle interplay between the critical singularities of the non-zero tem- 
perature transition and those of the T = quantum critical point. 

As in the Ising case, it is useful to take a Hamiltonian point of view. 
In this case (see Section 3.2 for justification) the associated model is the 
following Hamiltonian of O(N) quantum rotors on the sites i of a regular 
d dimensional lattice (we will consider general N > 2 for completeness): 



where the sum (ij) is over nearest neighbors, J > is an overall energy 
scale, and g > is a dimensionless coupling constant. The N component 
vectors Hi are of unit length, nf = 1, and represent the orientation of the 
rotors on the surface of a sphere in iV-dimensional rotor space; as we will 
see later, n — ► (j) a in the mapping to the CQFT S. The operators Liuv 
(fj, < v, fi, v = 1 . . .TV) are the iV(./V — l)/2 components of the angular 

momentum Li of the rotor: the first term in Hr is the kinetic energy of 
the rotor with 1/g the moment of inertia. The different components of fti 
constitute a complete set of commuting observables and the state of the 

system can be described by a wavefunction ^{fii). The action of Li on <I> 
is given by the usual differential form of the angular momentum 




(ij) 



(64) 




(65) 
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The commutation relations among the Li and Hi can now be easily de- 
duced. We emphasize the difference of the rotors from Heisenberg-Dirac 
quantum spins: the components of the latter at the same site do not com- 
mute, whereas the components of the rii do. 

While there are no direct realizations of quantum rotors in nature, they 
are a useful model of the low energy physics of a number of experimental 
systems. As we will see below, a quantum rotor should be considered as a 
model for the physics of a pair of antiferromagnetically coupled Heisenberg 

spins. The field n represents to the staggered antiferromagnetic order pa- 
's—* 

rameter, while L is the total magnetization. Thus quantum rotor models 
describe the spin fluctuations of the spin-ladder compounds [32], a rotor 
representing the spins on each rung of the ladder. For technical reasons 
not explored here, the rotor models also describe integer spin chains, and 
a large class of Heisenberg antiferromagnets in higher dimensions. Finally, 
rotor models are also useful in double-layer quantum Hall systems [33], 
where superexchange effects lead to an antiferromagnetic coupling between 
the electrons in the two layers [34]. These systems (excluding the double- 
layer quantum Hall system for which the reader is referred to the literature) 
will be described in more detail in the following subsections. 

There is a strong analogy between the rotor Hamiltonian Hr in (64) 
and the Ising Hamiltonian Hi in (5). We will be looking at the transition 
between a magnetically ordered state with (ft) ^ and O(N) symmetry 
broken, and a quantum paramagnet in which correlations of n are short 
ranged. As in the Ising model, it is the exchage term, proportional to J, that 
favors the ordered state, while the 'kinetic energy', proportional to Jg leads 
to fluctuations in the orientation of the order parameter and eventually to 
loss of long-range order. The similarity between the two models will also 
be apparent in the strong (large g) and weak coupling (small g) analyses 
in the following section. 

3.1. LIMITING CASES 

The pictures which emerge in the following two perturbative analyses are 
expected to hold on either side of a quantum critical point at g = g c . 

3.1.1. Strong coupling g 3> 1 

At g = oo, the exchange term in Hr can be neglected, and the Hamilto- 
nian decouples into independent sites, and can be diagonalized exactly. The 

<— > 

eigenstates on each site are the eigenstates of L 2 ; for N = 3 these are the 
states 



\£,m) i £ = 0,1,2,..., -£<m<£ 



(66) 
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and have eigenenergy Jg£(£ + l)/2. Compare this single site spectrum with 
that of a pair of Heisenberg-Dirac spin 5 quantum spins with an antiferro- 
magnetic exchange K; for a suitably chosen K we get the same sequence of 
levels and energies but with a maximum allowed value of £ = 25. Assuming 
this upper cutoff is not crucial for the low energy physics, we can use a 
single quantum rotor is an effective model for a pair of spins. 

The ground state of Hr in the large g limit consists of the quantum 
paramagnetic state with £ = on every site: 



Compare this with strong coupling ground state (7) of the Ising model. 
Indeed, the remainder of the strong coupling analysis of Section 2.1 can 
borrowed here for the rotor model, and we can therefore be quite brief. The 
lowest excited state is a 'particle' in which a single site has £ = 1, and this 
excitation hops from site to site. An important difference from the Ising 
model is that this particle is three-fold degenerate, corresponding to the 
three allowed values m. The dynamic susceptibility has a quasiparticle pole 
at the energy of this particle, and odd particle continua above the three 
particle threshold. 

3.1.2. Weak coupling, g <C 1 

At g = 0, the ground state breaks O(N) symmetry, and all the Hi vectors 
orient themselves in a common, but arbitrary direction. Excitations above 
this state consist of 'spin waves' which can have an arbitrarily low energy. 
This is a crucial difference from the Ising model, in which there was an 
energy gap above the ground state. The presence of gapless spin excitations 
is a direct consequence of the continuous O(N) symmetry of Hr: we can 
make very slow deformations in the orientation of (n), and get an orthogonal 
state whose energy is arbitrarily close to that of the ground state. Explicitly, 
for N = 3, and a ground state polarized along (1, 0, 0) we parametrize 



where 1 7ri | , | vr2 1 <C 1, and look at the linearized equations of motion for 7Ti, 
7T2. A standard calculation then gives harmonic spin waves with energy lo = 
ck (c is the spin wave velocity) ; their wavefunctions can then be constructed 
using harmonic oscillator states. 

3.2. CONTINUUM THEORY AND LARGE N LIMIT 

To obtain the path integral representation of the quantum mechanics of Hr, 
we interpret the Hi as the co-ordinates of particles constrained to move 



|0> = U\£ = 0,m = 0) i 



(67) 



n(x,t) = (l,7Ti(x,t),7r 2 (x,t)) 



(68) 
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on the surface of a sphere in A dimensions: we can then simple use the 
standard Feynman path integral representation of single-particle quantum 
mechanics. After taking the continuum limit, such a procedure gives 



Vn(x, r)5(n 2 (x, 



1) exp 



d d xC 



2cg 



dn\ 2 /dn\ 2 
dxiJ KdrJ 



(69) 



Here c ~ J a is a velocity which will turn out to be the spin wave velocity, 
and we have set h = ks = 1- The prefactor of A is for future convenience. 
The coupling constant g ~ ga d ~ x has the dimensions of (length)^ -1 ; we will 
not use the original g in Hr further in this discussion, and we will drop 
the tilde in g from now. The above action is valid only at long distances 
and times, so there is an implicit cutoff above momenta of order A ~ 1/a 
and frequencies of order cA. Our main interest here shall be the universal 
physics at scales much smaller than A. It is now also apparent that if we 
convert the fixed length field n to the "soft" spin cf) a , the action in (69) 
maps to the model S. 

A useful tool in the analysis of (69) is the limit of a large number of 
components of n i.e. expansion in 1/A [7]. For N > 3, the N = oo solu- 
tion already contains the proper dilineation of all the crossovers and phase 
transitions for all d. The description of the dynamic properties however re- 
quires a fairly subtle analysis of the fluctuations. As all of these issues have 
been discussed at length in the literature, here we will merely set-up the 
framework of the N = oo theory, and then proceed to a physical discussion 
of the properties of the various regimes for d = 1, 2, 3. 

The N = co solution is quite easy to set up, at least in the phase without 
long range order in the order parameter n; we will not discuss the ordered 
phase [7] here. We rescale the n field to 



n 



An, 



(70) 



and impose the h = N constraint by a Lagrange multiplier, A. The action 
is then quadratic in the n field, which can then be integrated out to yield 



J DA exp 



l/T 



< :2 i)j i) 2 ■ i\) ^ J ' dT J d d x\ 

V) 

The action now has a prefactor of A, and the A = oo limit of the functional 
integral is therefore given exactly by its saddle point value. We assume that 
the saddle-point value of A is space and time independent, and given by 
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iX = a 2 . The saddle-point equation determining the value of the parameter 
a 2 is 

a j^L tT i = 1 m 

(2vr) d ^ c 2 k 2 +to 2 +a 2 eg' { ' 

where the sum over ui n extends over the Matsubara frequencies u n = 2mrT, 
n integer. It is also not difficult to see that the retarded dynamic suscep- 
tibility x{k,w), which is the two-point correlator of the n field, is given 
by 

X( fc > W ) = - 2fc2 f Cg { N -X\2 i 2' ( 73 ) 

crk z — (oj + iby + a z 

at N = oo, with 5 a positive infinitesimal. The Eqns (72,73) are the central 
results of the N = oo theory; in spite of their extremely simple structure, 
these equations contain a great deal of information, and it takes a rather 
subtle and careful analysis to extract the universal information contained 
in them [7]. We will not go into these technical details here, but will be 
satisfied by describing the different physical regimes predicted by (72,73) 
and other analyses of Hr. 



3.3. DYNAMICS IN ONE DIMENSION: APPLICATION TO SPIN-GAP 
COMPOUNDS 

Solution of (72) in d = 1 shows that there is an energy gap, A, for all positive 
values of g. Even for very small values of g, where one might expect long- 
range order, and the weak coupling solution of (3.1.2) to apply, quantum 
fluctuations are strong enough to always drive the system paramagnetic. 
For small g, the energy gap behaves like 

A ~ cAe" 2 ^ 9 . (74) 



This behavior can be confirmed by a more sophisticated renormalization 
group analysis [35]. The implication of this important result is that the 
large g analysis of Section 3.1.1 actually applies for all values of g in d = 1. 

The finite T crossovers can also be deduced by a solution of (72) [36], 
and are sketched in Fig 9. The quantum critical point of Fig 3 has effec- 
tively been moved to g = 0, and there are now only T distinct limiting 
physical regimes: (i) T <C A, where strong quantum fluctuations create a 
singlet paramagnetic ground state, and we have a few thermally excited el- 
ementary excitations — for N = 3 these are a triplet of spin 1 particles with 
energy (A 2 +c 2 k 2 ) 1 / 2 at momentum k, and correspond (crudely) to break- 
ing a singlet valence bond between two neighboring spins in the underlying 
antiferromagnet; (ii) A <C T <C cA, the high T regime of the continuum 
theory, where quantum fluctuations are marginally subdominant to thermal 
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Figure 9. Crossovers of Hr (Eq. (64)) in d — 1 for N > 3. The phase boundary is at 
A ~ T, with A given by (74). 



fluctuations [37] (by a factor of l/ln(T/A)), and we can locally describe 
the system in terms of a doublet of left/right circularly polarized spin- 
waves about an ordered state — however, interactions of thermally excited 
spin waves lead again to a paramagnetic state. 

Below we will present an exact theory for dynamic and transport prop- 
erties in the T <C A regime [38]. There is as yet no complete theory for 
dynamics in the T> A regime, and its description remains an important 
open problem. 

The basic approach followed for T <C A is very similar to that taken 
for the quantum-disordered region of the Ising chain in Section 2.3.2. The 
central difference here is that the quasi-particle excitations have an addi- 
tional spin label, m = 1 . . . N . We will present our discussion below for the 
physical case N = 3, although the generalization to arbitrary ./V is immedi- 
ate. This label on the quasiparticles is associated with the conserved O(N) 
charge, which is the total spin J2i ^i- We shall restrict our attention here to 
the dynamics associated with the transport of this conserved charge: this 
will be described by obtaining the exact long-time form of the correlator 

C af3 (x,t) = (L a (x, t)L p (0,0)) (75) 

For the case N = 3, the second-rank antisymmetric tensor indices \iv on L 
have been replaced by a single vector index a or (3. The method described 
below can also be used to obtain dynamic correlators of the n field: this 
will be discussed elsewhere [39]. 

There are two key observations that allow our exact computation for 
T -C A. The first, as in the Ising chain, is that as the density of particles 
~ e -A/T^ anc j ^heir mean spacing is much larger than their thermal de- 
Broglie wavelength ~ c/(AT) -1 / 2 ; as a result the particles can be treated 
semiclassically. The density of particles with each spin m (m = —1,0, 1 for 
N = 3) is now given by the expression (39), and the total density therefore 
equals 3p. The second observation is that collisions between these particles 
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Figure 10. Two particle collision described by the 5- matrix (76). The momenta before 
and after the collision are the same, so the figure also represents the spacetime trajectories 
of the particles. 

are described by their known two-particle 5-matrix [40] , and only a simple 
limit of this 5-matrix is needed in the low T limit. The r.m.s. velocity 
of a thermally excited particle vt = c(T / A) 1 / 2 , and hence its 'rapidity' 
~ vt/c <C 1. In this limit, the 5-matrix for the process in Fig 10 is [40] 

5 m?,m| = (-^^rmm'^m^- (76) 

In other words, the excitations behave like impenetrable particles which 
preserve their spin in a collision. As in the Ising chain, energy and mo- 
mentum conservation in d = 1 require that these particles simply exchange 
momenta across a collision (Fig 10). The (—1) factor in (76) can be in- 
terpreted as the phase-shift of repulsive scattering between slowly moving 
bosons in d = 1. Indeed, the simple form of (76) is due to the slow motion of 
the particles, and is not a special feature of relativistic continuum theory: 
it can be shown that (76) also holds for lattice Heisenberg spin chains in 
the limit of vanishing velocities. 

As in Sections 2.3.1 and 2.3.2, we represent C(x,t) as a 'double time' 
path integral, and in the classical limit, stationary phase is achieved when 
the trajectories are time-reversed pairs of classical paths (Fig 11). Each 
trajectory has a spin label, m, which obeys (76) at each collision. The 
label, m, is assigned randomly at some initial time with equal probability, 
but then evolves in time as discussed above (Fig 11). We label the particles 
consecutively from left to right by an integer k; then their spins are 
independent of t, and we denote their trajectories Xk{t). The longitudinal 
correlation C33 is then given by the correlators of the classical observable 

L 3 (x,t) = ^2m k 5(x - x k (t)) (77) 
k 

in the classical ensemble defined in Section 2.3.1 and above. Now because 
the spin and spatial co-ordinates are independently distributed, the cor- 
relators of rrik and x^ factorize. The are uncorrelated, and therefore 
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1 ►x 

Figure 11. A typical set of particle trajectories contributing to C(x,t). Each trajectory 
represents paths moving both forward and backward in time, and the ( — 1) phase at each 
collision is neutralized by its time-reversed contribution. The particle co-ordinates are 
Xk{i), with the labels k chosen so that Xk(t) < xi(t) for all t and k < I. Shown on the 
trajectories are the values of the particle spins nik which are independent of t in the low 
T limit. 



(rrikme) = (2/3)5ke We therefore have 

C 33 (x - x ', t -t') = ^Yl ( 5 ( x ~ x k{t))5{x' - x k (t'))) (78) 
6 k 

This result involves the self two-point correlation of a given particle k, 
which follows a complicated trajectory in the way we have labeled the 
particles (e.g. the trajectory of the —1 in Fig 11). Fortunately, precisely this 
correlator was considered three decades ago by Jepsen [41] and a little later 
by others [42]; they showed that, at sufficiently long times, this correlator 
has a Brownian motion form. Inserting their results into (78), we obtain 
the results presented below. 

An important property of the results is that they can written in a 're- 
duced' scaling form [7], much like that found in Sections 2.3.1 and 2.3.2. 
The characteristic spatial length (L x ) and time (Lt) can be chosen to be 

1 1 / A \ 1//2 

L * = r P L < = ^U?r) ■ (79) 

Notice L T ~ ce A / T /\/AT is the mean spacing between the particles , and 
Lt ~ e^l T /T is a typical time between particle collisions. Our final result 
is 

C aP (x,t) = 6p 2 F (M, B.) 5 a/3 (80) 

where F is a universal scaling function. A complete expression for F is given 
in Ref [38]. Here we note the important limiting forms. For short times F 



34 



has the ballistic form 

F(x,t) ^ e~^ p /ty/n, (81) 

which is the auto-correlator of a classical ideal gas in d = 1, and holds for 
|t] <C |x| <C 1. In contrast, for \i\ » 1, \x\ it crosses over to the diffusive 
form 

In the original dimensionful units, (79) and (82) imply a spin diffusion 
constant, D s , given exactly by 

c 2 e A / T 

D, - ifj-. (83) 

To the best of our knowledge, this is first exact result for the spin diffusivity 
in any paramagnetic spin system. 

The above long-time diffusive form of the magnetization correlations 
have important consequences for NMR experiments on one-dimensional 
compounds with a spin gap. The comparison with experiments requires 
generalization of the above theory to include an external magnetic field H. 
This has been discussed elsewhere [38]. From this computation one finds 
that the longitudinal nuclear spin relaxation rate, 1 jT\ , is given in the weak 
field limit H < A, L^ 1 by 

_1_ FT X _ TAe- 3A / 2r f^T _ e' A / T 

where T is a nuclear hyperfine coupling, and Xu is the paramagnetic spin 
susceptibility (in the rotor model, this is the linear response to a field 
H under which Hr — > Hr — H^ i Li^). For experimental comparisons, a 
striking property of the above, pointed out to us by M. Takigawa, is that the 
low T activation gaps for \/T\ (Ajm) and % (A x ) satisfy Aw^/A^ = 3/2. 

Values for the activation gaps A]y Tl and A x have been quoted for a 
number of experimental systems [43, 44, 32], and it has consistently been 
found that A-im ls larger than A x . The above theoretical results provide 
a natural explanation for this puzzling experimental trend. This is further 
explored in Fig 12, which was kindly provided to us by H. Yasuoka. 

We have also compared [39] the detailed H and T dependence of the 
1/Ti measurements on Ref [43] with the predictions of our theory. These 
measurements have enough dynamic range to explore the ballistic-diffusive 
crossover described above, and the results are in good agreement with the- 
oretical predictions. 
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Figure 12. Measured values of the activation gaps A 1; / Tl and A x for a number of spin 
gap systems. The one-dimensional compounds are enclosed by full rectangles, while others 
have dashed rectangles (these are either two-dimensional or consist of isolated dimers). 
The full line is the theoretically predicted ratio of 3/2 and is in good agreement with the 
trends in the d — 1 systems. The other systems have a ratio close to unity, and this is 
expected to be the case when elastic impurity-scattering of the quasiparticles dominates. 

3.4. TWO DIMENSIONS 

Solution of the N = oo equation (72) in d = 2 yields a phase diagram 
essentially identical to that shown in Fig 3 for the d = 1 quantum Ising 
model. A renormalization group analysis [21] shows that these large N 
results are qualitatively correct for the physical case of iV = 3. A detailed 
description of the physical properties in the different regions of the phase 
diagram of Fig 3 for the d = 2 0(3) rotor model has been obtained in the 
1/N expansion [7]. We will not describe this analysis here as the physical 
properties are quite similar to those for the d = 1 Ising model, which we 
have discussed in some detail in Section 2. Further, the results of the large 
N theory have been reviewed elsewhere [3, 45]. 

We do note one important difference from the d = 1 Ising model, how- 
ever. The ordered phase (g < g c ) of the d = 2 0(3) rotor model has gapless 
spin-wave excitations above the ground state, unlike the Ising model which 
had a gap. Consequently, we need a different energy scale to characterize the 
deviation of the ground state from the quantum critical point g = g c . A con- 
venient choice is the spin stiffness, p s , which has the dimensions of energy in 
d = 2. The non-zero T properties of the rotor model will now be universal 
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functions of the dimensionless ratio p s /T for g < g c . The renormalized- 
classical region is therefore present for T < p s , while the "High T" region 
is T > p s . There is an energy gap, A, towards a triplet of quasiparticle 
excitations in the quantum-disordered phase present for g > g c , and the 
nonzero T crossovers in this case are very similar to those of the d = 1 
Ising model: the physics is "quantum-disordered" for T < A, and "High 
T" for T > A. 

The d = 2 0(3) rotor model is expected to describe the low energy 
properties two dimensional quantum antiferromagnets whose low energy 
states are well described by fluctuations of a collinear antiferromagnetic 
spin ordering. The 5 = 1/2 spin- 1/2 square lattice antiferromagnet, real- 
ized in La2Cu04 and related insulators, is an important example of such a 
system, and we now discuss the relationship of its properties to those of the 
rotor model. By a comparison with low temperature behavior of the corre- 
lation length observed in neutron scattering measurements, Chakravarty et 
al. [21] convincingly argued that the ground state of this model has long- 
range Neel order, i.e. it maps onto the rotor model with a coupling g < g c . 
Their analysis was restricted to the renormalized classical region of Fig 3 
i.e. T <C p s , where static properties reduce to those of an effective clas- 
sical rotor model. A theory for the dynamic properties in this region was 
presented in Ref [46] by postulating classical equations of motion for the 
classical rotor model. 

Next, we turn to a question [47] which arises naturally from the struc- 
ture of Fig 3: does the 5 = 1/2 square lattice antiferromagnet crossover 
into the "High T" regime at temperatures T > p s ? The existence of such 
a crossover is subject to the condition that it occur at temperatures small 
enough so that the mapping to rotor model remains valid. Strong evidence 
for the "High T" regime was found in Ref [47] by a comparison of the 
1/N expansion of the rotor model with numerical and experimental mea- 
surements of the uniform susceptibitlity. The magnitude and temperature 
dependence of the 1/T\ [48] and I/T2G [49] nuclear magnetic relaxation rate 
measurements on La2Cu04 were also found to be in good agreement with 
theory [7, 50]. Convincing evidence for the renormalized-classical to "High 
T" crossover was presented by Sokol et al. [51] and Elstner et al. [52] using 
a high temperature series expansion for a number of static correlators of the 
5=1/2 square lattice antiferromagnet. (Elstner et al. [52] also examined 
the 5=1 square lattice antiferromagnet, and found no evidence for the 
"High T" regime; this indicates, as expected, that the effective value of g 
for 5 = 1 is significantly smaller than g c .) It was also argued early on [47] 
that the T dependence of the correlation length, which was central to the 
conclusions of Chakravarty et al. , would not exhibit a clear signature of 
the "High T" regime. This point was subsequently reiterated by Greven et 
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al. [53]. 

A more complete understanding of the physics of the "High T" regime 
should emerge from a careful study of the wavevector, frequency, and tem- 
perature dependence of the imaginary part of the staggered dynamic spin 
susceptiblity. This should obey a scaling form [6] much like that obtained 
earlier for the Ising model in (61): 

///i x A - ( hck Tluj \ 

w,» ) = —* x {—,—y (84) 

where A is a non-universal amplitude and 77 is the anomalous dimension 
of ft field at the quantum-critical point (77 ~ 0.03 for d = 2, N = 3). 
Detailed theoretical results for the form of the universal scaling function <3? x 
have been obtained [6, 7] in the 1/N expansion. Unfortunately there are no 
existing results for x"(k, u>) in La2Cu04 at T > p s , as the neutron scattering 
intensity becomes quite small at higher temperatures. However there are 
other two-dimensional square lattice antiferromagnets with a smaller p s [54] 
which are good candidates for such a study. 

This is a good point to make some speculative remarks about the doped 
antiferromagnet L^-xSr^CuO,^ The 1/Ti measurements of Ref [48] have 
the striking property that while the low T values of 1/Ti are quite strongly 
doping dependent, at higher T 1/Ti approaches a value which is both tem- 
perature and doping-concentration independent. This appears to receive 
a natural explanation within the framework of a rotor model quantum- 
critical point. Assuming the main effect of doping is merely to shift the 
value of the effective coupling constant g, the properties of the "High T" 
region of Fig 3 are determined primarily by the value of T alone, and are 
insensitive to the precise value of g: indeed, in this region, the system does 
not 'know' whether g is smaller or larger than g c . Further, the T depen- 
dence is predicted to be 1/Ti ~ T T ' [47, 7], and given the small value of 77, 
this is essentially T- independent. Lastly, the coefficient of T v has also been 
estimated [7] and is consistent with the experimental value. 

We turn now to other antiferromagnets in d = 2 which have a collinear 
Neel ordered ground state. An interesting system in which it is possible 
to tune the value of the effective coupling g across g c is the double layer 
square lattice S =1/2 antiferromagnet [55]. This model consists of spin- 1/2 
Heisenberg spins on two adjacent square lattices, with an intralayer anti- 
ferromagnetic exchange J and an interlayer antiferromagnetic exchange K. 
The ratio Kj J acts much like the dimensionless coupling g, with the large 
Kj J phase a gapped quantum paramagnet of singlet pairs of spins in oppo- 
site layers, and the small Kj J phase Neel ordered. Numerical simulations 
have been carried out on this model by Sandvik and collaborators [56], and 
the critical point g = g c identified rather reasonably well. A number of 
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Figure 13. Phase diagram for the 0(N) rotor model with for N > 2 and d = 3. A 
similar diagram applies for the cases with N = 1, 2 and d = 2,3. The dashed lines are 
crossovers, while the full line is the locus of finite temperature phase transitions. The 
shaded region is where the reduced classical scaling functions apply. 



universal amplitude ratios have been studied at this point, and all results 
are now in good agreement with the 1/N expansion on the 0(N) quan- 
tum rotor model. High temperature [52, 57] and strong coupling [58] series 
expansions on the double-layer model also reach a similar conclusion. 

An exciting recent development has been the remarkably high precision 
study of quantum-critical properties in the 1/5-depleted square lattice by 
Troyer et al. [59]. This is another model in which g can be tuned across g c 
by varying the ratio of exchange constants. This particular lattice has the 
additional feature that the mapping to the rotor model is expected to be 
incomplete: there are additional Berry phases generated by the dynamics 
of the Heisenberg spins undergoing a "hedgehog" tunneling event [60]. It 
was argued [7] that such perturbations should be "dangerously irrelevant" , 
and that they do not modify any of the universal scaling predictions of the 
rotor model. The results of Troyer et al. are in excellent agreement with 
the critical properties of the rotor model [59], thus supporting the neglect 
of the Berry phases. 



3.5. THREE DIMENSIONS 

Solution of the N = oo equation (72) now yields the phase digram shown in 
Fig 13 [45]. The new feature is presence of a finite temperature phase tran- 
sition: this lies entirely within the low T region for g < g c and is not a direct 
property of the quantum critical point. The well-studied classical critical 
behavior appears only within the shaded region of Fig 13. The quantum- 
critical scaling functions contain this classical criticality as "reduced scaling 
functions" in much the same way that in the d = 1 Ising model the classical 
function (46) was contained in the more general quantum function in (35). 
This phenomenon can be studied with the 1/N expansion [45], and has 
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also recently been described in an expansion in e = 3 — d [8]. Also, because 
d = 3 is the upper-critical dimension of S, there are logarithmic violations 
of scaling near the quantum-critical point. 

Normand and Rice [61] have proposed an interesting recent experi- 
mental realization of such a three-dimensional quantum critical point in 
LaCu02.5- This is a spin-ladder compound in which the ladders are mod- 
erately strongly coupled in three dimensions. By varying the ratio of the 
intra-ladder to inter-ladder exchange it is possible to drive such an anti- 
ferromagnet across a d = 3 quantum critical point separating Neel ordered 
and quantum paramagnetic phases. The uniform susceptibility has a T 2 
dependence at intermediate T, which is characteristic of the "High T" re- 
gion in d = 3. The entire T dependence of the uniform susceptibility has 
been compared with large ./V and quantum Monte Carlo simulations of the 
quantum-critical model [59], with good agreement. 

4. Quantum relaxational transport in two dimensions 

Transport properties of the model S were considered earlier in Section 3.3. 
There we determined the spin diffusion constant in the quantum-disordered 
regime (T <C A) of the N = 3 model in d = 1. Here, following Ref [63], we 
shall describe transport for N > 2 in d = 2 in the "High T" region (Fig 1 
and 3), where the dynamics is quantum-relaxational (Fig 8) . The methods 
discussed here can also be applied to the quantum-disordered region of the 
N > 2, d = 2 models, but we shall not present that extension here. Let us 
also recall that there is one region for which there is no quantitative theory 
of transport phenomena: the "High T" region for N > 3 in d = 1; this 
region should be experimentally accessible in S = 2 spin chains. 

We shall present our discussion using the language of the N = 2 system, 
although closely related results apply to all N > 2. The N = 2 model 
describes a superfluid-Mott insulator transition in a lattice boson model 
with short-range repulsive interactions [64]. The complex superfluid order 
parameter VP is related to the field 4> a by 

* = (0i + ifa)/V2 (85) 

(and similarly to n). It therefore serves a starting point for understand- 
ing superfluid-insulator transitions in disordered thin films [65, 66] and 
Josephson junction arrays [67] and the transitions among quantum Hall 
states [68, 69, 70]. 

Rather than describing the transport in terms of the diffusivity, it is 
more convenient here, both theoretically and experimentally, to characterize 
it using the conductivity, a; the two quantities are, of course, connected by 
the Einstein relation. We consider the response of S to a spatially varying 
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"chemical potential", fx, (this is the field that couples to the conserved 
charge for N = 2; for the N = 3 case we called it the "magnetic field") 
under which the rotor model Hamiltonian transforms as 



Here Q is the charge of the field. Note that for N = 2 there is only one 
generator L which therefore does not carry any vector indices. In an open 
system, the expectation value of the current (this is the current associated 
with the conserved charge) will be proportional to the gradient of and 
the proportionality constant is a. We will be interested here in the uj and T 
dependence of a at the critical point g = g c , as that controls the value of o 
in the "High T" region. We begin by writing down the expected scaling form 
this must satisfy. For this we need the scaling and engineering dimensions 
of a. Rather general arguments [71, 72, 73, 68] show that in d = 2, a has 
scaling dimension zero, and engineering dimensions of the quantum unit of 
conductance Q 2 jfi. This leads to the scaling form 



where S is a universal scaling function whose properties we will describe 
here. 

It is clear from (87) that the experimentally measured d.c. conductivity 
is determined by the universal number X(0). In contrast, any computation 
of the conductivity performed at T = 0, in which the uj — > limit is 
taken subsequently will determine the number S(oo). It was argued quite 
generally by Cha et al. [72] and Wallin et al. [74] that £ is in independent 
of u/T, which therefore also implies that E(0) = S(oo). As a result, there 
were a number of analytic computations [71, 72, 75, 76, 77, 78, 79, 80, 81] 
and an exact diagonalization study [82] of the value of S(oo) in a variety of 
models which display a quantum-critical point in d = 2. Quantum Monte 
Carlo studies [72, 74, 83, 84, 85] determined the conductivity by analytic 
continuation from imaginary time data at uj = 2irnTi, with n > 1 integer; 
as the smallest value of uj is 2nTi, these studies were effectively also in the 
regime |a;| 3> T. 

The picture of quantum relaxational dynamics we have reviewed in this 
paper makes it quite clear that the underlying assumption of these works is 
incorrect: £ is in fact, not independent of u/T. As we have discussed here, 
the characteristic property of the "High T" region is that there is a phase 
relaxation time t v ~ fi/ksT. Dynamic order parameter fluctuations also 
carry charge, and therefore inelastic collisions between thermally excited 
charge-carrying excitations will lead to a transport relaxation time r tr . As 




(86) 



z 
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Figure 14- A sketch of the expected form of the real part, E', of the universal scaling 
function E appearing in the scaling form (87) for the conductivity, as a function of 
ZJ — hLo/ksT. There is a Drude-like peak from the inelastic scatterings between thermally 
excited carriers at 57 of order unity. At larger W, there is a crossover to the collisionless 
regime where E' —* E(oo) . 



the typical energy exchanged in a collision is ksT, r tr is also of order t v , 
and therefore 

v ~ ¥■ < 88 > 

r tr h 

The missing coefficient in (88) is a universal number whose value will de- 
pend upon the precise definition of l/rt r . It is perhaps worth noting ex- 
plicitly here that (88) holds for all values of the dynamic exponent z. The 
relationship (88) could be violated if interactions were dangerously irrele- 
vant at the quantum critical point, a possibility we shall not discuss here. 
Now general considerations [86] suggest that there are two qualitatively 
different regimes of charge transport at non-zero frequencies 

(I) uJTtr <C 1 (hu <C ksT): the hydrodynamic, collision-dominated, inco- 
herent regime, where charge transport is controlled by repeated, inelastic 
scatterings between pre-existing thermally, excited carriers; the conductiv- 
ity should exhibit a 'Drude' peak as a function of frequency. 

(II) ujT tT 3> 1 (frw 3> ksT): the high frequency, collisionless, phase-coherent 
regime, where the excitations created by the external perturbation are solely 
responsible for transport, and collisions with thermally excited carriers can 
be neglected. 

These physical arguments lead to a suggested form for the function £(uJ) 
shown in Fig 14. We have assumed that S(0) > S(oo), and this will be 
the case in the specific calculation for S discussed below. The minimum 
at co ~ ksT jU also appears in the calculation below, but it is not clear 
whether that is a general feature. 

It is also interesting to consider the implication of Fig 14 for a(uj,T — > 
0). This is illustrated in Fig 15 which plots the form of (h/Q 2 )a(uj,T — > 0) 
in d = 2: its value at ui = is given by S(0), while for all uj > it equals 
E(oo). Note the difference from Fermi liquid theory, where the Drude peak 
becomes a delta function with non-zero weight as T — > 0. In the present 
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Figure 15. Universal form of the conductivity a(co, T — » 0) in d = 2; the vertical scale is 
measured in units of h/Q 2 . Only the w = value is given by the universal number E(0). 
For all w > 0, {h/Q 2 )a = E(oo). 



case, the weight in the Drude-like peak vanishes like ~ T as T — > 0, and 
reduces to the single point cj = where the conductivity is given by S(0). 

We note that the above discussion implies results that are rather re- 
markable from a traditional quantum transport point of view. There have 
been a number of previous situations in which charge transport properties 
have been found to be universally related to the quantum unit of conduc- 
tance, e 2 /h; these include the quantized Landauer conductance of ballistic 
transport in one-dimensional wires, and the universal conductance fluc- 
tuations of mesoscopic metals [87]. However in all previous cases, these 
universal properties have arisen in a "phase-coherent" regime, i.e. they are 
associated with physics at scales shorter than the mean distance between 
inelastic scattering events between the carriers. For the case of a d = 2 
quantum critical point discussed above, the universal number £(oo) is as- 
sociated with quantum coherent transport, and is therefore the analog of 
these earlier results. In contrast, the value of S(0) is controlled by repeated 
inelastic scattering events, and therefore the d.c. transport is clearly in 
what would traditionally be identified as the "incoherent" regime. Never- 
theless, we have argued above that S(0) is a universal number, and the d.c. 
conductance remains universally related to (Q = e) 2 /h. 

Let us now turn to the computation of £ for the specific model S. The 
T > transport properties of S in d = 2 were considered briefly by Cha 
et al. [72]. They concluded that the continuum model S had a = oo for all 
T > 0, and that additional lattice umklapp scattering effects were needed 
to degrade the current carried by thermally excited carriers, and to yield a 
a < oo. However, in Section 3.3 we a found finite diffusivity D s for N = 3, 
d = 1 at T > 0, and we claim that a similar result holds here for N = 2, 
d = 2. The basic point can be made quite simply by a glance at Fig 11. Note 
that in this d = 1 model, the transport of momentum is indeed ballistic, 
and is represented by the straight line trajectories which extend for all time. 
However, when we consider the transport of charge, we have to follow the 
motion of the iV-valued O(N) index carried by each trajectory: this motion 
is seen in Fig 11 to be rather complicated, and was shown in Section 3.3 to 
be diffusive. There is thus an essential difference between momentum and 
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charge transport when the excitations carry a O(N) flavor index, and it 
appears to have been overlooked in Ref [72]. 

We will outline the recent computation of £ in an expansion in e = 
3 — d [63]. For small e the quartic interaction, uq, in S reaches a universal 
(to leading order in e) fixed point value [1] 

24vr 2 , , 

n = —e (89) 

in the vicinity of the quantum critical point. So the scattering of the el- 
ementary excitations, and the non-ballistic motion of the charge current, 
only occurs with an amplitude of order e. Such a weak-coupling 5-matrix 
should be compared with the complementary strong-coupling S'-matrix in 
(76) for d = 1: in the latter case we had a maximal scattering situation 
as momentum and spin labels interchanged their mutual pairing in every 
collision. Also in contrast to the real-space formulation of transport taken 
for d = 1 in Section 3.3, the weak-coupling situation here is best analyzed 
in momentum space. We wish to work with same charge eigenstates we 
considered in Section 3.3, and so we express the field * in the associated 
normal modes 



* (X '* } = / 72^72T (MMK kx W_(M)e-* kx ) , (90) 



where Sk = {k 2 + m 2 ) 1 / 2 , with m a renormalized T-dependent mass [63]. 
The operators a± annihilate particles with 0(2) charges ±1 (in Section 3.3, 
where we considered 0(3), there were three such particles). Next, transport 
theory requires the momentum space distribution functions 

/ A (M) = (4(M)a A (k,t)), (91) 



with A = ±1, in terms of which a quantum Boltzmann equation for trans- 
port in a weakly-scattering system can be obtained [86]. We will not display 
the explicit form of this equation, or discuss its solution here: the reader 
should consult Ref [63]. That analysis leads to an explicit solution for £ 
which has the qualitative form of Fig 14. We quote quantitative results for 
the two limiting values [63, 78]: 

E(0) _ ™»+... 

ol-2d_l-d/2 

s '<^°°> = ^fW2H 1+0<£) )" (92) 

to leading order in e. Notice the singular e dependence of the first result: 
this is signal of a "boundary layer" in S of width uj ~ e 2 T which constitutes 
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the "Drude peak" of Fig 14. This boundary layer explicitly demonstrates 
a claim made in Section 1: the operations of expansion in e and of ana- 
lytic continuation in to do not commute. The structure of the higher-order 
corrections to the first result in (92) is quite complex, and was generally 
discussed in Ref [8]. For the superfluid insulator transition in d = 2 with 
Q = 2e, (92) implies a universal d.c. conductivity 

4e 2 4e 2 
<j(0) = 2^S(0)— » 1-037— (93) 

This result is remarkably close to the value, 4e 2 //i (argued [75, 76] to be 
related to self-duality), and to the results of many experiments [66]. Indeed, 
it was conjectured in Ref [63] that 27rS(0) = 1 exactly for N = 2 in d = 2. 
Definitively establishing this self-duality however requires techniques other 
than expansion in e = 3 — d, or 1/N, as it is only possible precisely at d = 2, 
N = 2. Finding such techniques is an important challenge for future work. 

5. Conclusion 

We have reviewed the universal, non-zero temperature dynamic properties 
of the CQFT S for a number of representative cases, and discussed their 
experimental implications. Among the cases omitted, there are some that 
show phenomena qualitatively different from those discussed here: 
(i) The case N = 2, d = 1 has a line of critical points and a Kosterlitz- 
Thouless end point at T = 0. Its paramagnetic phase is similar to the 
N = 3, d = 1 model discussed in Section 3.3. 

(h) For d > 3 there is violation of hyperscaling, which leads to explicit 
dependence of finite T behavior on the bare value of uq and qualitatively 
different crossovers [8]. 
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